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We analyze theoretically the electronic properties of Aharonov-Bohm rings made of graphene. We 
show that the combined effect of the ring confinement and applied magnetic flux offers a controllable 
way to lift the orbital degeneracy originating from the two valleys, even in the absence of intervalley 
scattering. The phenomenon has observable consequences on the persistent current circulating 
around the closed graphene ring, as well as on the ring conductance. We explicitly conflrm this 
prediction analytically for a circular ring with a smooth boundary modelled by a space-dependent 
mass term in the Dirac equation. This model describes rings with zero or weak intervalley scattering 
so that the valley isospin is a good quantum number. The tunable breaking of the valley degeneracy 
by the flux allows for the controlled manipulation of valley isospins. We compare our analytical 
model to another type of ring with strong intervalley scattering. For the latter case, we study a ring 
of hexagonal form with lattice-terminated zigzag edges numerically. We find for the hexagonal ring 
that the orbital degeneracy can still be controlled via the flux, similar to the ring with the mass 
confinement. 

PACS numbers: 73.23.-b, 73.23.Hk, 73.23.Ra, 81.05.Uw 



I. INTRODUCTION 

Graphene offers the remarkable possibility to probe 
predictions of quantum field theory in condensed mat- 
ter systems, as its low-energy spectrum is described by 
the Dirac- Weyl Hamiltonian of massless fermionsi How- 
ever, in graphene, Dirac electrons occur in two degenerate 
families, corresponding to the presence of two different 
valleys in the band structure - a phenomenon known as 
"fermion doubling" . This valley degeneracy makes it dif- 
ficult to observe the intrinsic physics of a single valley in 
experiments, because in many cases the contribution of 
one valley to a measurable quantity is exactly cancelled 
by the contribution of the second valley. A prominent 
example that single-valley physics is interesting, is the 
production of a fictitious magnetic field in a single val- 
ley by a lattice defect or distortion.^ The field has the 
opposite sign in the other valley, so its effect is hidden 
when both valleys are equally populated. Another exam- 
ple is the existence of weak antilocalization in diffusive 
graphene, which is destroyed by intervalley scattering,'^ 
Therefore, from a fundamental point of view, it is de- 
sirable to find a feasible and controlled way to lift the 
valley degeneracy in graphene. From a more practical 
point of view, the lifting of the orbital degeneracy is es- 
sential for spin-based quantum computing in graphene 
quantum dots^, which is a promising direction of future 
research because of the superior spin coherence proper- 
ties expected in carbon structures. 

Here, we show that the confinement of electrons in 
graphene in an Aharonov-Bohm (AB) ring (see Fig. 1) 
provides a conceptually simple way to achieve a con- 
trolled lifting of the valley degeneracy. We find that the 
ring confinement - which breaks the effective time rever- 
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FIG. 1: A circular graphene ring of radius a and width W 
subjected to a magnetic fiux $ threading the ring. 

sal symmetry (TRS) in a single valley in the absence of 
intervalley scattering^ - leads generically to a lifting of 
the valley degeneracy controllable by magnetic flux. To 
demonstrate this, we choose an analytical model with a 
smooth ring boundary described by a mass term in Sec- 
tion II. Such a mass term might be generated in a real 
system by the influence of the lattice of the substrate 
on the band structure of the sample, see Refs. Hl^ for 
concrete examples of such an effect. However, our con- 
clusions are not restricted to the mass confinement but 
potentially hold for any boundary that preserves the val- 
ley isospin as implied by general symmetry arguments)^ 
We show that the signature of the broken valley degen- 
eracy is clearly visible in the persistent current and the 
conductance through the ring. It is further illustrated 
how to use the lifting of the valley degeneracy with flux 
to manipulate and measure the valley isospin. 

In Section III, we compare our analytical model for a 
smooth boundary to a system where intervalley scatter- 
ing is strong. For this purpose, we study a ring of hexago- 
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nal shape with zigzag edges where intervalley scattering is 
induced at the corners of the hexagon. We calculate the 
spectrum numerically in a tight-binding approach and 
find that the orbital degeneracy can still be tuned by the 
magnetic flux, similar to the analytical model. We test 
this ability against a small distortion of the 6-fold sym- 
metry of the ring and find small avoided crossings at zero 
flux. 



II. RING WITH SMOOTH BOUNDARY 

In this section, we analyze in detail the spectral prop- 
erties of a graphene ring subjected to a magnetic flux 
and its signatures in persistent current and conductance 
through the ring assuming a smooth confinement induced 
by a space-dependent mass term in the Dirac equation. 
We also discuss how to address the valley degree of free- 
dom in such a ring structure. 
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A. Spectrum 



The graphene ring with valley degree of freedom t = ± 
is modeled by the Hamiltonian (?i = c = 1) 

Hr = Ho + TV{r)a,, (1) 

where we use the valley isotropic form^ for Hq — v{p + 
eA) • a with p — —id/dr, —e being the electron charge, 
V the Fermi velocity and (7i=x,y,z are the Pauli matri- 
ces. The vector potential is A = ($/27rr)e^ with $ the 
magnetic flux threading the ring, see Fig. 1. The term 
proportional to az in Eq. ([1]) is a mass term confining 
the Dirac electrons on the ring. Introducing polar coor- 
dinates, the Hamiltonian i/o is written as 

Ho^r^if) = —iv{cosipax + smipay)dr 

If ^\ 

— iv{cos (fiCTy sin ipax)— [dip + i— j . (2) 

The angular orbital momentum in the z-direction is Iz = 
—id^ and $o = 27r/e. The two valleys r = ± decouple 
and we can solve the spectrum for each valley separately, 
H-rtpT = Eipr- We note the similarity of Ht to a ring 
with the Rashba interactionjiSiiii However, an important 
difference is that for graphene, the confinement potential 
acts on the (pseudo)spin, whereas for Rashba interac- 
tion the confinement potential is spin- independent. As 
a consequence, the confining potential in Eq. ([T]) couples 
the pseudo-spin components and breaks effective TRS 
{p —p, a — > — (t) even in the absence of a flux <I>i^ 



Since Hr commutes with 
■0T are eigenstates of Jz, 



h + ^fz, its eigenspinors 



_ gi(m-l/2)ip 



(3) 



FIG. 2: Energy spectrum i5om {E > 0) as a function of mag- 
netic flux $ for various total angular-momentum values m 
(blue m > 0, red m < 0) using Eq. ^ with a/IF=10: (a) 
shows a single valley only r = -1-1. The dotted lines are the 
exact numerical evaluation of E using Eq. ([6]). In (b) we show 
the full spectrum including the other valley r = — 1 (dashed 
lines) . The restoration of ±$ symmetry in the combined spec- 
trum of both valleys and the lifting of the valley degeneracy 
at finite $ are clearly visible. 



with eigenvalues m, where m is a half-odd integer, 
TO — ±i,±|,... The radial component Xt{t) = 
(xT(»^),X2(^)) satisfies Hr{r)xT{r) = Exr{r) with 



HT-{r) — —ivcTxdr +TV{r)az + 



I f m — 

V(T„- I ^ 







(4) 

where we have defined m = m -I- ($/$o)- 

For V{r) — 0, xK'^) and xJl*") are solutions to Bessel's 
differential equation of order to — ^ and rn + ^, respec- 
tively. Therefore, the eigenspinor for Hr(r) with energy 
E and V{r) — can be written as 



Xt 



-br 



m — 4- ^' ' 



(5) 

where Hi ' (p) are Hankel functions of the (first, sec- 
ond) kind and the dimensionless radial coordinate is 
p = \E\r/v. The coefficients a^. and h^- are determined 
by the boundary condition of the ring induced by V{r) 
(with V{r) — > -|-oo outside the graphene ring). We use 
the infinite mass boundary condition ij^T = T(n^ • (T)'ipr 
where n._L = (— sin ip, cos ip) a,t r = a + ^ and with op- 



posite sign at r 



W 5J. 
2 ^ 



Here a is the ring radius and 
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W its width (see Fig. 1). Eliminating the coefRcients Or 
and br gives the energy eigenvalue equation z — z* with 



(1) 



(1) 



(6) 



which is equivalent to cj) — im with cj) the phase of 
z and n an integer. In Eq. ([5]), we have abbreviated 
pi = \E\{a - \)lv and p2 = \E\{a + \)lv. To obtain 
an analytical approximation of the spectrum, we use the 
asymptotic form of the Hankel functions for large p, in- 
cluding corrections up to order 1/ p^}'^ This indeed is the 
desired hmit as p = \E\r/v ~ \E\a/v oc a/W » 1 when 
the ring radius is much larger than its widtbi^ and leads 
to the following energy eigenvalue equation 



\E\ = ^. 
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(7) 



An iteration of Eq. ([7]) by replacing |£'| on the right-hand- 
side of the equation by the first (leading) term of \E\ gives 
the energy eigenvalues (neglecting terms of 0[(Vl^/a)^]) 



e: 



(n+1) 



(8) 



In Eq. 



, En — v{n + ^)n/W, n = 0,1,2,..., and 

A„ = {v I of /2£n- These energy eigenvalues are plot- 
ted as a function of flux for n = and different values of 
m (its half-odd integer values reflect the 7r-Berry phase 
of closed loops in graphene) in Fig. 2. Fig. 2(a) shows the 
energy levels for one valley, r = -1-1. It is clearly visible 
that E'^{rn) ^ E'^{—rn), since effective TRS is broken by 
the confinement. In Fig. 2(b), the spectrum of both val- 
leys is shown with full lines for r = + 1 , and dashed lines 
for r = -1. At $ = 0, E'^{m) = E-'^{-m) as it should 
be, since real TRS is present at zero magnetic field. Cru- 
cially, however, at finite $, E~^ ^ E~ in general, showing 
that the valley degeneracy is indeed lifted since effective 
and real TRS are broken. If n ^ 1, the term oc rm 
in Eq. ^ becomes suppressed and the valley degener- 
acy is restored, correctly predicting that the spectrum 
is insensitive to the boundary condition if 27r/q„ <C W, 
where g„ = 7r(n -|- \)/W is the transverse wave number. 
We show next that a broken valley degeneracy results 
in observable features in the persistent current and the 
conductance through the ring. 

B. Persistent current and conductance 

The persistent current in the closed ring is given at zero 
temperature by j = — Tlnm '^^nml'^^ where the sum 
runs over all occupied states. In Fig. 3, we show the per- 
sistent current as a function of number of electrons on 
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FIG. 3: Persistent current as a function particle number 
(including spin) and flux $ for n = [E > Qi): (a) in- 
cludes only valley t = whereas (b) includes both valleys. 
Curves for difi'erent A*' are displaced with dashed horizontal 
lines defining j = for each curve. The broken valley degen- 
eracy is clearly visible in (b) via two substructures of length 
A<I> — 2$o/7r and A$ = (1 — (2/7r))$o whereas (a) predicts 
a non-zero persistent current at $ = due to effective TRS 
breaking in a single valley. 



the ring N (including spin) and magnetic flux relative to 
the half-filled band. (We subtract the contribution to the 
persistent current that arises from all states with E <Q.) 
The persistent current is periodic in $ with period $0- 
In Fig. 3(a), only one valley, r = +1 is considered and a 
finite persistent current at <I> = is predicted. Therefore, 
a non-zero persistent current at zero fiux detects valley 
polarization. In Fig. 3(b), we show the case of equal pop- 
ulation of both valleys. Then, the persistent current as a 
function of flux is zero at <I> = 0, but shows a substruc- 
ture (kinks at $ 0) within one period directly related 
to the broken valley degeneracy at finite flux. We note 
that this substructure is due to the linear term in m of 
the spectrum Eq. ([5]) which is prominent within the first 
few transverse modes n which can host many electrons 
N. 

In Fig. 4, we plot the conductance through the ring 
weakly coupled to leads as a function of Fermi energy Ep 
(or gate voltage) assuming a constant interaction modeli^ 
with charging energy U At <& = 0, the conductance ex- 
hibits a four-fold symmetry due to spin and valley degen- 
eracy. A finite flux breaks the valley degeneracy which is 
observable via a splitting of the conductance peaks mov- 
ing with magnetic flux, see Fig. 4. 
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FIG. 4: Ring conductance assuming a constant interaction 
model with charging energy U for the first 12 electrons in 
the conduction band (-E > 0): At $ = (dashed), the con- 
ductance shows a four-fold symmetry as a function of Fermi 
energy Ep in the leads due to spin- and valley-degeneracy. 
At finite magnetic flux {$/^o = 0.1, full line), the conduc- 
tance peaks shift due to breaking of the valley degeneracy. 
Each peak is labelled with the filling factor N at this specific 
resonance. 



C. Valley qubit 

We now turn to the question of how to make use of 
the broken valley degeneracy in order to directly address 
the valley degree of freedom in graphene experimentally 
(valleytronicslS) . The valley degree of freedom forms 
(in principle) a two-level system that can be represented 
by an isospin |-|-) for valley t = +1 and |— ) for valley 
T = — 1. We point out that the graphene ring weakly 
coupled to current leads could be used to investigate the 
relaxation and coherence of such valley isospins. Close 
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FIG. 5: Valley qubit. The lowest four levels (for n = and 
a/VK=10) are shown as a function of flux. Blue and red 
(dashed) lines correspond to valleys r = -1-1 and r = — 1, 
respectively. The dotted lines take into account small level 
mixing leading to anticrossings. The flux <I> is used to switch 
from |-|-) to the crossing point with the new eigenstates being 
superpositions of |-|-) and |— ) as indicated in the figure. 



to a degeneracy point of two levels belonging to different 
valleys (e.g. at $ = 0), Eq. ([S]) predicts a valley splitting 
of states with fixed m-values controllable by flux, simi- 
lar to the Zeeman-splitting for electron spins in a mag- 
netic field. In semiconductor quantum dots, such pairs of 
spin-split states can be addressed via electron tunneling 
from/to leads weakly coupled to the quantum dot and 
can be used for read-out of single spinsii or measuring 
their relaxation (Ti) timeJ^ The graphene ring could be 
used in very much the same way to measure the intrinsic 
valley isospin relaxation time Ti in graphene as well as 
the valley isospin polarization. 

In Fig. 5, we show the situation when some small level- 
mixing leads to avoided crossings of valley-split states 
near the degeneracy point $ = 0. Such valley mixing nat- 
urally appears through boundary roughness of the ring 
or atomic defects in the bulk. Using the magnetic flux 
as a knob, we can sweep the system from a |-t-) ground- 
state level, filled with one electron, to a superposition 
(1+) + l~))/v^ and further to a |— ) groundstate, see 
Fig. 5. Such a situation can be used to produce Rabi 
oscillations of the valley isospin states by tuning the sys- 
tem fast (non-adiabatically) from |-|-) to the degeneracy 
point where the spin will oscillate between |-|-) and |— ) in 
time: cos(A<)|-|-) — ism{At)\—), where 2A is the energy 
splitting at the degeneracy pointji^ 

We expect that this qubit is rather robust if intervalley 
scattering is weak, since time-reversal symmetry assures 
the (approximate) degeneracy of states from different val- 
leys at zero flux Indeed, the Hamiltonian Eq. (1) has 
the same spectrum in both valleys at zero flux, indepen- 
dent on the shape of the mass potential V {x , y) .'^'^ This 
means that we do not rely on a special symmetry of the 
confining potential (like the circle discussed here). In 
addition, long-range disorder will also not lift the valley 
degeneracy. 



D. Valley isospin-orbit coupling 



In an open ring geometry with adiabatic contacts to 
leads, new interesting coherent rotations of the valley 
isospin occur while propagating along the ring. The lin- 
ear term in m in Eq. ^ can be thought of as a valley 
isospin-orbit coupling term, since the valley isospin t cou- 
ples to the orbital motion m. A general incoming spinor 
is a superposition of spinors belonging to different val- 
leys. Due to the valley isospin-orbit coupling, the angular 
momentum m [determined by the incoming (continuous) 
energy E and the applied magnetic flux $ via Eq. ([S])] will 
be different for the two valleys. Consequently, the spinor 
in Eq. ([3]) will pick up different phases expiirmp) for the 
two valleys while propagating along the ring thereby ro- 
tating the valley isospin in a transport experiment. 



5 



III. SPECTRUM FOR A HEXAGONAL RING 
WITH ZIGZAG EDGES 

Here, we compare our analytical model with the in- 
finite mass boundary described in Section II, to a ring 
with strong intervalley scattering. We numerically in- 
vestigate the spectrum of a ring of hexagonal form with 
zigzag edges as shown in Fig. 6. (Electrical conduction 
through this geometry was studied in Ref. [2l|.) In a 




FIG. 6: A hexagonal ring with zigzag edges of inner radius 
Ta and outer radius and flux through the hole "I>. Here, 
Ta = 3\/3' and rt 



Q\fdil with I the lattice spacing 



zigzag nanoribbon, the valley isospin is a good quantum 
number, i.e. the zigzag boundary does not mix valleysi^ 
Since two neighboring arms of the ring are rotated by 
60° with respect to each other, the roles of the A and 
B sublattices are interchanged in subsequent segments. 
Explicitly, this means, that if a zigzag edge is terminated 
on a A side, it will be terminated on a B side at a neigh- 
boring arm of the ring. Equivalently, in the reciprocal 
(k-) space, this implies that equivalent states of subse- 
quent zigzag nanoribbon segments are lying in opposite 
valleys. This necessarily induces intervalley mixing at 
the corners between two subsequent zigzag nanoribbon 
segments. This mixing is very strong in the lowest mode 
of the ring, where the direction of motion and the valley 
is tightly coupled in each arm of the hexagonal ringi^i^ 
(the zigzag edge is therefore another example where ef- 
fective TRS in a single valley is broken). An electron 
wave, approaching a corner of the hexagon in one valley, 
is either transmitted into the next arm, or reflected back 
into the same arm. In both cases, the valley index has to 
flip. 

We investigate the spectrum of such a ring numerically 
in a tight-binding approach with Hamiltonian 



(9) 



The hopping element in the presence of a magnetic flux is 
tij — — t exp[— i(27r/<I'o) /r""' dr ■ A] where A is the vector 
potential and ci = are the on-site energies. The vector 
potential is chosen as A = {Ax^ 0, 0) with 
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FIG. 7: Plot (a) shows an energy band of the hexagonal ring 
(see Fig. 6) with ring dimensions Ta = 7^3/ and n = 14^3; m 
the lowest mode < -E ^ 0.34f.™ The levels are grouped into 
bands containing 6 levels. The top most and the lowest level 
in each band is non-degenerate whereas the middle four levels 
are two-fold degenerate at zero flux. This degeneracy is lifted 
by the flux through the ring. In (b) we contrast the perfect 
crossing of two levels at zero flux (red dots) with anticrossed 
levels (blue circles) induced by the addition of one unit cell 
to each of two parallel arms of the ring [we have shifted the 
energy axis for the asymmetric case (blue circles) by -|-3- 10~^t 
for better comparison]. 



where ya{x) = min[r'a, V3(Lc/2 - |a;|)], with Lc = 
'^fa/Vi for the perfect hexagon (shown in Fig. 6) and 
Lc — Ara/Vi + I for a hexagon ring with one unit cell 
added to the top and bottom arm. This represents a 
uniform magnetic field B inside the ring hole and zero 
outside. The spectrum as a function of magnetic flux 
$ = 2y/?>r^B is shown in Fig. 7 in an energy win- 
dow which lies well within the lowest mode of a zigzag 
nanoribbon of width VF = rf, — ra.~ A band of levels in 
the lowest mode is shown in Fig. 7(a). Within that mode, 
the spectrum follows a clear pattern which is observed for 
generic values of and r^. It consists of bands separated 
by energy gaps. Each band hosts six levels. The top and 
bottom level is non-degenerate with dE/d^ = at zero 
flux (corresponding to standing waves). The other four 
levels are two-fold degenerate at $ = with a broken de- 
generacy at finite flux. These levels correspond to right 
and left-going states in the ring. 
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We remark that this level pattern reflects the scatter- 
ing off a periodic array of six scatterers subjected to pe- 
riodic boundary conditions.^'* It is to be noted that the 
orbital degeneracy of the hexagonal ring can be tuned by 
the flux, similar to the ring with the smooth conflncment 
discussed in Section II. If the 6-fold rotational symmetry 
of the ring is broken, the crossings at zero flux become 
slightly avoided as is shown in Fig. 7(b) (blue circles) 
where we have added one unit cell to two of the paral- 
lel arms of the ring (this corresponds to a length change 
of the arms by about 5%). This shows that our results 
are also relevant for rings with a sligthly reduced sym- 
metry. Note that the sensitivity of the level crossing at 
zero flux to the ring geometry is consistent with strong 
intcrvalley scattering where time-reversal symmetry does 
not protect the degeneracy at zero flux. 

IV. CONCLUSION 

We have analyzed the Aharonov-Bohm effect in 
graphene rings. We have investigated two different ring 
systems - a ring with a smooth boundary (with zero or 
weak intervalley scattering) and a hexagonal ring with 
zigzag edges. For the ring with a smooth boundary, 
the combined effect of the effective time reversal symme- 
try (TRS) breaking within a single valley induced by a 
smooth boundary and the applied magnetic flux - break- 
ing the real TRS - gives us at hand a controllable tool 
to break the valley degeneracy in such rings. We have 



shown that this effect of a broken valley degeneracy by 
flux is revealed in the persistent current and in the ring 
conductance. This tool could be useful for spin-based or 
valley-based quantum computing. The presence of a de- 
generate pair of levels from different valleys at zero flux 
is assured by time-reversal symmetry in the absence of 
intervalley scattering. Therefore, the proposed effect is 
not sensitive to the actual geometry of the ring. 

We have also considered the opposite case of strong in- 
tervalley scattering by investigating numerically a hexag- 
onal ring with lattice-terminated zigzag edges. Here, 
strong intervalley scattering is induced by the corners 
of the ring at low energies. We found that the orbital de- 
generacy of graphene can still be tuned by the flux similar 
to the ring with a smooth boundary. This effect, how- 
ever, relies on a certain degree of symmetry of the ring 
as we show by sligthly distorting the ring. We therefore 
conclude that the orbital degeneracy in graphene rings 
can be controlled with an Aharonov-Bohm flux in rings 
with zero or weak intervalley scattering and in systems 
with strong intervalley scattering if the ring possesses an 
(approximate) geometric symmetry. 
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